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ABSTRACT 

Spiral galaxies that move through the intracluster medium lose a substantial amount 
of their gas discs due to ram pressu re stripping. The recent observations of NGC 4388 
bv lOosterloo fc van Gorkori] l)2005() reveal a tail of stripped gas of ~ 100 kpc behind 
the source galaxy. We present first 3D hydrodynamical simulations of the evolution of 
such ram pressure stripped tails. 

We find that if the ICM wind does not vary significantly over a period of a few 
100 Myr, subsonic galaxies produce a tail with regular features similar to a von-Karman 
vortex street. In this case, the tail widens systematically by about 45 kpc per 100 kpc 
distance behind the source galaxy. The widening rate is independent of the galaxy's 
inclination for a large range of inclinations. For supersonic galaxies, the tail is more 
irregular than for subsonic ones. 

The tail observed for NGC 4388 is narrower than the tails in our simulations. 
Reasons for this difference may be additional physical processes such as heat conduc- 
tion or viscosity. In addition, we conclude that the observed S-shape of this tail is not 
due to von Karman oscillations, because this galaxy is likely to move supersonically. 
A reason for the observed shape may be motions in the ambient ICM. 

Finally, we discuss implications for the distribution of metals in the ICM due to 
ram pressure stipping. 



Key words: 

medium 



galaxies: spiral ~ galaxies: evolution ~ galaxies: ISM ~ intergalactic 



1 INTRODUCTION 

Ram pressure stripping, i.e. the removal of a galaxy's gas 
disc due to its motion through the intracluster medium 
(ICM), has been stu died mainly with r egard to the effect 
on the galaxv [e.g. lAbadi et al.l fl999l: lOuilis et al.l l2000l:_ 
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ISchiilz fc Struckll200ll: IVollmer et alJl200ll: iMarcolini et alJ 
l20n,'j:lR,oediger fc Hensle jl200.4 iRocdigcr fc Bruggen "200^. 
It has been shown that the interaction between the ICM 
and the galactic interstellar medium (ISM) can remove 
a significant amount of gas from galaxies near cluster 
centres. Thus, ram pressure stripping is a good candi- 
date to explain observations s uch as the HI deficiency of 
cluster disc galaxi es (see e.g. ICavatte et al.l 1199(1 . Il994l : 
ISolanes et 3111200 J) and the truncated star forming dis cs 
in t he Virgo cluster feoopmann fc Kennevlll99^ . l2004al lhl: 
iKoo pmaim ct al. 2005). A distinct feature of ram pressure 
stripped galaxies is that their gas discs are distorted or trun- 
cated while their stellar discs appear undisturbed, because 
ram pressure stripping affects only the geiseous components 
of the galaxy. Such examples have been observed, e.g. NGC 
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4522 ("Kennev fc KooDmann"l999','200lVKonncv et al."2004 
FVoUmer et al...2004) NGC 4548 (.VoUmer ct al. 1999) , NGC 
4848 ijVollmer et al.ll2001^ . In addition to the effect on the 
galaxy, ram pressure stripping plays an important role in 
the evolution of the ICM because th e gas lost by galaxies 
enriches the IC M with metals (see e.g. lSchindler et alJl2005| : 
iDomainko et aL.2005. and references therein). 

A further example of ram pressure stripped galaxies is 
the Virgo spiral NGC 4388. It has a small distance to the 
cluster centre and a high radial ve locity (^ 1400kms~^) 
with respect to the cluster mean (see lVollmer fc Huchtmeien 
and references therein) . There is strong evidence that 
it has sufi'ered ram pressure stripping in i ts recent his- 
tory. It was observed to be HI deficient by ICavatte et alJ 
l|l99nh . Moreover, both, ionised and neutral gas have been 
found at a few kpc distance (also sec Vollmor fc Huchtmeiral 
12003 for further references). VoUmcr & Huchtmcicr ( 20031) 
reported the detection of atomic gas at a projected dis- 
tance of more than 20 kpc from the galactic centre. This 
extraplanar gas is interpreted as the tail of gas stripped 
from NGC 4388. On the basis of sticky particle simulations, 
IVollmer fc Huc htmeier ( 2003) concluded that the main mor- 
phological features and velocity structure of NGC 4388 and 
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the extragalactic gas clouds can be explained by ram pres- 
sure stripping. 

Recently, lOosterloo fc van GorkomI ||2QQ^ presented 
observations of a ~ 120 kpc long tail of HI gas associated 
with NGC 4388, and also suggest that this tail is due to 
ram pressure stripping of this galaxy, either in the ICM of 
the Virgo cluster or in the halo of the nearby elliptical galaxy 
M86. Two prominent features of the tail are its flaring width 
and its slight S-shape. 

Although ram pressure stripping has been studied by 
several groups, only very few investigations of th e re- 
sulting gas tai l s are available. IStevens et alJ ^1999^ and 
lAcreman et al.l (|200.'j ) simulated ram pressure stripping of 
hot gas halos from spherical galaxies. They predicted that 
bow shocks and a tail of hot material will be observable 
in X-rays. Both groups used cylindrical 2D hydro codes for 
their simulations. 

To our knowledge, so far no simulations of tails of cool 
gas stripped from disc galaxies have been presented. How- 
ever, the physics that determines the dynamics of such tails 
is complex and interesting: 

• According to basic hydrodynamics, any object that 
moves through a fluid produces a wake. The characteris- 
tics of this wake depend on th e Reynolds numb er (see any 
hydrodynamics text book, e.g. iBatcheloJ l200Cfl . For small 
Reynolds numbers (Re < 50), the wake is thin. For moder- 
ate Reynolds numbers (Re ~ 100), features similar to the 
von Karman vortex street develop. The wake's width can 
reach several times the width of the body. For high Reynolds 
numbers (Re ~ 1000), the wake becomes turbulent. Espe- 
cially for moderate and high Reynolds numbers, the wake 
widens signiflcantly with increasing distance to the object. 
The value of the Reynol ds number in the ICM is a matter of 
ongoing debate (see e.g. iRevnolds et al.ll2005ll . with a large 
range of values suggested. 

• The shght S-shapo of the tail of NGC 4388 present in 
the observation of Oostcrlo o fc van G orkom (2005) is rem- 
iniscent of a von Karman vortex street. However, classical 
von Karman vortex streets are usually generated behind "2D 
bodies", like cylinders, that force the flow into a transla- 
tional symmetry along the axis of the body. In contrast, the 
flow past a galaxy is a full 3D flow and thus possesses a more 
complex wake structure. 

• The stripping process alone may lead to a widening of 
the galaxy's tail because the ICM is forced to flow around 
the galactic gas disc. Therefore, the galactic gas can aquire 
velocity components perpendicular to the main ICM wind 
direction. 

• Additions to the dynamics in the galactic wake/tail 
could arise from the rotation of the disc gas. As the stripped 
galactic gas moves away from the galaxy, the radial gravita- 
tional acceleration ceases to operate on the stripped gas. In 
the extreme case of a face-on galaxy where the radial accel- 
eration is simply switched off, a stripped gas package would 
move linearly with v = Urot perpendicular to iTicM from that 
moment on. The ICM wind can accelerate the stripped gas 
package to maximal iiicm in wind direction. The superpo- 
sition of a motion with vicm in wind direction and Vrot in 
any direction perpendicular to this would lead to a widen- 
ing angle of a = 2 arctan(?;rot/i'iCM). This process should be 
inclination-dependent. Of course, stripped gas has to pass a 



transition region where the radial gravitational acceleration 
ceases to operate. 

• The inclination angle between the galactic rotation axis 
and the galaxy's direction of motion may influence the wake 
due to geometrical reasons alone, i.e. the effective cross- 
section of the "obstacle" in the ICM wind changes. Further- 
more, the flow patterns of subsonic and supersonic cases are 
expected to differ because in supersonic cases a bow shock 
develops. 

Most of the aspects listed above require a 3D treatment 
of the problem. Here, we present flrst 3D hydrodynamical 
adaptive-mesh simulations of the evolution of tails of ram 
pressure stripped galaxies. 



2 METHOD 

We study the motion of the galaxy through the cluster in 
the rest frame of the galaxy. Therefore, its motion translates 
into an ICM wind flowing past the galaxy. 



2.1 Code 

The simulations we re performed with the FLASH code 
l|Frvxell et al.l l200nh . a multidimensional adaptive mesh 
refinement code. It solves the Riemann problem on a 
Cartesian grid using the Piecewise-Parabolic Method 
(PPM), which is particularly well suited to resolve shocks. 
The gas obeys a polytropic equation of state with an 
adiabatic index of 7 = 5/3. All boundaries but the inflow 
boundary are open. The simulations presented here are 
performed in 3D. In order to follow the evolution of the 
wake, we use a large simulation box of size (a;min,Kmax) x 

(ymin.ymax) X (Zmin , ^niax) = (- 162 kpC, 162 kpc) X 

(-64 kpc, 260 kpc) x (-162 kpc, 162 kpc) for subsonic 
cases and ( — 121.5 kpc, 121.5 kpc) x (—64 kpc, 260 kpc) x 
(—121.5 kpc, 121.5 kpc) for supersonic ones. The galactic 
centre is located at (a;gai, j/gai, Zgai ) = (0,0,0). The ICM 
wind is flowing along the y-axis into the positive direction. 
We use 6 levels of reflnement, which leads to an effective 
resolution of 316.5 pc and an effective number of grid cells 
of (1024)'' for the subsonic cases. However, at a distance of 
more than 64 kpc behind the galaxy (behind=downstream) 
we limit the effective resolution to 633 pc (5 levels of 
refinement) to limit computational costs. For the same 
reason, the number of refinement levels is restricted to 5 for 
the supersonic runs, where a multitude of shocks and waves 
causes massive refinement (see Fig. I^J. The refinement 
criteria are the standard density and pressure gradient 
criteria. The galactic disc is always refined. 

In order to be able to identify the galactic gas after it is 
stripped from the galaxy, we use a "dyeing" technique: The 
FLASH code offers the opportunity to advect mass scalars 
along with the density. We utilise one of these mass scalars, 
/, to contain the fraction of galactic gas in each cell. Initially 
this array has the value 1 in the region of the galactic disc 
and elsewhere. As a result, every cell where / > contains 
a certain amount of gas that has originally been inside the 
galaxy. At each timestep, the quantity fp gives the local 
density of galactic gas. 
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Table 1. Galaxy model parameters. 
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2.2 Model galaxy 

We model a massive spiral galaxy with a flat rotation curve 
at 200 kms"^. It consists of a gas disc, a stellar disc, a stellar 
bulge and a dark matter (DM) halo. For the gravitational 
potential of the stellar disc, bulge and DM halo, we use the 
following analytical descriptions: 

Stellar disc: Plummer-Ku zmin disc, se e 

iMivamoto fc Nagail (jlflTfJ^ or iBinnev fc Tremainel ljl987l) . 
Such discs are characterised by their mass, M*, and radial 
and vertical scale lengths, a* and 6,, respectivel y. 

Stel lar bulge: Spherical Hernquist bulge (see iHernaui'sl 
Il99,^) . In case of a spherical bulge, the gravitational poten- 
tial, "1>, depends on radius, r, as <I>(r) = £i^^£buigc ^ where 
Afbuige is the mass of the bulge, rbuige the scale radius and 
r the spherical radius. 

DM halo: The spherical model of lBurker'3 (|l^9.'^^ . including 
the self-scaling relations, i.e. the DM halo is characterised 
by the radial scale length, tdm, alone. For the equation of 



the analytical potential see also iMori fc BurkertI ^ 

The self-gravity of the gas is neglected as the gas contributes 
only a small fraction to the overall galactic mass. The pa- 
rameters of our galaxy model are summarised in Table Q 
In order to prevent steep density gradients in the galactic 
plane and in the galactic centre, the gas disc is described by 
a softened exponential disc : 



p{R,Z) = 



27r(Igas ^gas 



0.5^ sech ( — ) sech ( 7^ ) . (1) 



The coordinates {R, Z) are galactic cylindrical coordinates. 
The radial and vertical scale lengths are flgas and bgas, 
respectively. For R > Ogas and \Z\ > fegas, this density 
distribution converges towards the usual exponential disc 
p{R,Z) = 2na-'''l — exp(-i?/agas) exp(-|Z|/&gas). For the 

gas 

corresponding exponential disc, Mgas is the total gas mass. 
We chose Mgas such that in the outer regions the gas disc 
converges to an exponential gas disc with 10% of the stel- 
lar disc mass, i.e. Mgas = O.IM*. The integrated gas mass 
amounts to 6 • 10® Mq. Given the density distribution in the 
disc, its pressure and temperature distribution are set such 
that hydrostatic equilibrium with the ICM is maintained in 
the direction perpendicular to the disc plane. In radial direc- 
tion, the disc's rotation velocity is set so that the centrifugal 
force balances the gravitational force and pressure gradi- 
ents. We have cut the gas disc smoothly to a finite radius of 
26 kpc by multiplying the density distribution p(R, Z) with 
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Figure 1. Radial profiles of the density, p, pressure, p, and ro- 
tation velocity, frot, in the galactic plane for the initial model. 
Also the radial profile for the projected ISM surface density, S, 
is shown. 



Table 2. ICM wind parameters. Ram pressure pram, ICM 
density picMi ICM wind velocity "UicMi inclination angle i. 
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0.5[1 + cos{tt{R - 20kpc)/6kpc)] for 20 kpc < i? ^ 26 kpc. 
Figure Q shows radial profiles of density, surface density, 
pressure and rotation velocity for the initial model. 

We measure the inclination angle, i, of the galaxy as the 
angle between its rotation axis and the ICM wind direction. 
Thus 0° corresponds to to a face-on motion of the galaxy, 
and 90° to edge-on. 



2.3 ICM conditions 

For a first experiment, we expose model galaxies to con- 
stant ICM winds. We list Pram, vicm, Picm and inclination 
angle, i, in Table |5| We chose the ICM temperature such 
that the sound speed is 1000 km s-^. We refer to the two 
ram pressures as medium (pram = 6.4 • 10-^2 grgcm-'^) and 
strong (piam ~ 6.4 • IQ-^^ ergcm-^) ram pressure. For both 
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ram pressures, we test a subsonic (uicm = 800kms~^, Mach 
number 0.8) and a supersonic (iiicm = 2530 km s~^, Mach 
number 2.53) case. 

For the subsonic runs, we start with the ICM at rest and 
then increase the inflow velocity during the first 50 Myr. For 
the supersonic runs, this procedure takes too long until a 
true supersonic flow including a bow shock develops. There- 
fore, for supersonic runs, we set vicm everywhere in the ICM 
already initially. In this procedure, the simulations starts a 
bit violently, but after a bit less than 100 Myr a true bow 
shock has developed. This difference in initialisation leads 
to a temporal offset between subsonic and supersonic runs. 
A time t in a supersonic run corresponds to approximately 
t + 50 Myr in a subsonic run. 



3 RESULTS 

3.1 Snapshots in slices 

Figure|21shows the gas density distribution in a slice through 
the simulation box parallel to the y-z-plane at a; = for dif- 
ferent timesteps for the case of medium ram-pressure and 
supersonic ICM flow. In addition to the bow shock, the 
stripped material induces a rich shock structure in the re- 
gion behind the bow shock. The stripped material fragments 
strongly due to turbulence and Rayleigh- Taylor-instability. 

3.2 Projected gas densities 

In this section, we show projected galactic gas densities 
along two lines-of-sight for our simulation runs at different 
timesteps (Figs. |3| to |7|l. The galactic gas is identified by 
the mass scalar described in Sec. 12.11 In all figures in this 
section, the left column shows projections along the a;-axis, 
while the right column shows projections along the z-axis. 
These are the two directions perpendicular to the direction 
of motion (which is along the y-axis). The time is denoted 
in the upper left panel of each snapshot. Each figure cap- 
tion denotes the ram pressure, Mach number and inclination 
angle of the run. 

All projections of the medium ram-pressure cases show 
gas tails with column number densities around 10^^ cm~^. 
The tails of the high ram-pressure cases have a column den- 
sity roughly a factor of 10 higher. Typical tail widths are 
between 60 and 100 kpc. 

In the cases with medium ram-pressure, the galaxy can 
retain a tail for several 100 Myr and the tail can also evolve 
during this time. In the cases with strong ram-pressure, the 
complete gas disc is stripped after a few 100 Myr, so that 
one cannot speak of a galactic tail any more. Instead, the 
bulk of stripped gas moves downstream as a stretched gas 
cloud which is few 100 kpc long. The fact that in the high 
ram-pressure cases much more gas is lost in a shorter time 
leads to a distribution of more mass into a smaller volume 
and thus higher column densities in the tails. 

3.3 Velocity structure in the tail 

Figures |H1 to 1121 show the velocities of the stripped gas in 
the wake. For a random, ISM density weighted subset of 
grid cells, the velocity components in y-, x- and z-direction 
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Figure 2. Slice through simulation box in the y-z-p\ane for the 
supersonic run with medium ram- pressure and i = 30°. The local 
gas density is shown colour-coded for different time steps. The 
refinement is demonstrated in the first plot, one square shows the 
size of one block= 8^ grid cells. 



are plotted as a function of distance to the galaxy. In our 
simulations, Vy is the velocity component in wind direction, 
whereas and are perpendicular to the wind direction. 
The caption of each figure denotes the ram pressure, Mach 
number and inclination angle of the run. 

The component Vy reflects the acceleration of the 
stripped gas in the wind direction. The components perpen- 
dicular to the wind direction cannot result solely from galac- 
tic rotation. In the panels for the Vx- and u^-components, we 
have marked the amplitude of the galactic rotation velocity 
Vrot = ±200 kms~^. If the stripped gas were pushed straight 
into the wind direction, no velocities higher than this should 
occur in and Vz . But in all cases velocities higher than the 
galactic rotation occur. Especially in the 75° case, the galac- 
tic rotation induces only a small u^-component. Almost all 
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Figure 3. Projected galactic gas density for medium ram- 
pressure, Mach number 0.8, i = 30°. Lowest contour is 
10~®g cm~'^=6.0 ■ 10^® cm~^. The contour spacing is half an or- 
der of magnitude. The left column is for projection along the 
X-axis, the right column for projection along the 2-axis. The time 
of each panel is denoted in its upper left corner. The two numbers 
in the upper right corner denote the projection direction with re- 
spect to the galaxy: Oax is the angle between the line-of-sight 
and the galactic rotation axis. If the line-of-sight is projected into 
the grid's x-z-plane, the resulting line has an angle of ciproj with 
the negative part of the grid's z-axis. aax and cVproj E're given in 
degree. 



of the observed Vz must be due to the ICM flow. The excess 
in Vx and arises because the ICM is flowing smoothly 
around the galaxy. Thus, already the ICM flow near the 
galaxy contains velocity components in x- and ^-direction. 
In the supersonic case, the oblique bow shock also leads to 
non-zero Vx- and «z-components. The strongest increase of 
Vx and Vz is observed in the supersonic case. Additional ex- 
cess in Vx and is introduced by turbulence. 

The velocity width in the tail also differs in subsonic and 
supersonic cases. The overall width in the subsonic cases is ~ 
400-600 km s~^, whereas in the supersonic cases it is nearly 
twice as much. The velocity width is roughly independent 
of ram pressure. 

Interestingly, in the subsonic cases, oscillations appear 
in Vx and especially in Vz- Their origin is discussed below. 



Figure 4. Projected galactic gas density for medium ram- 
pressure, Mach number 0.8, i = 75°. Otherwise same as Fig. l3l 
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Figure 5. Projected galactic gas density for medium ram- 
pressure, Mach number 2.53, i = 30°. Otherwise same as Fig. l3l 
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Figure 6. Projected galactic gas density for strong ram-pressure, 
Mach number 0.8, i = 30°. Otherwise same as Fig. [3 In the lower 
panels, the bulk of stripped gas can be seen to move away from 
the galaxy (the galactic centre is located at (x, y, z) = (0, 0, 0) for 
all timesteps). 
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Figure 7. Projected galactic gas density for strong ram-pressure, 
Mach number 2.53, i = 30°. Otherwise same as Fig. ISl 



3.4 Tail structure 

The shape of the tail depends mainly on Mach number. Only 
in the subsonic cases regular structures similar to a von Kar- 
man vortex street can be produced (see Figs.j^and^J- The 
formation of a long regular tail also takes several 100 Myr. 
Therefore, less regular structure is produced in the strong 
ram-pressure run because the gas disc is already stripped 
completely after a few 100 Myr. Moreover, in Fig. El an os- 
cillating tail appears in the projection along the x-axis, but 
not along the 2;-axis. If we neglect the disc rotation for the 
moment, then the galaxy appears asymmetrical due to its 
inclination only in z-direction, but not in a;-direction. There- 
fore, the "von Karman oscillations" are more prominent in 
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Figure 8. Velocity components in the wake, for medium ram- 
pressure, Mach number 0.8, i = 30°. The top panels shows the 
velocity along wind direction (i/-axis), the middle and bottom 
panels show the velocity components in the two directions per- 
pendicular to the wind direction (x- and z-direction). Each col- 
umn is for one timestep. The time is denoted in the top panel of 
each column. The dashed vertical lines in the middle and bottom 
panels denote the amplitude of galactic rotation. 




10D 150 7QQ 750 



5P 1 go 150 7QQ 7S0 
y/kpc 



Figure 9. Velocity components in the wake, for medium ram- 
pressure, Mach number 0.8, i = 75°. Otherwise same as Fig. ISl 



^-direction, and hence can be seen best if the line-of-sight 
is the X-axis. Also the oscillations observed in the velocity 
components perpendicular to the wind direction Vx and Vz 
(see Sect. are connected to the "von Karman oscilla- 
tions" . 

The supersonic cases show an irregular tail structure. 
3.5 Flaring of the wake 

In addition to the von-Karman oscillations, in the subsonic 
cases the tail widens systematically although this effect is 
less prominent in the high ram-pressure case because the 
gas disc is already stripped completely at early times. The 
velocity components perpendicular to the wind direction are 
similar for, both, the medium and high ram-pressure case. 
Therefore, the tail is slightly narrower in the stronger ram- 
pressure case because here the stripped gas is accelerated 
faster. 

In the subsonic, medium ram-pressure cases, the widen- 
ing proceeds roughly linearly on the length scale of our sim- 
ulation box. The tail width ranges from a few lOkpc near 
the galaxy to 150 kpc at the end of the simulation box (dis- 
tance of 260 kpc behind the galaxy). The diagonal lines in 



© 2006 RAS, MNRAS 000.IT1ITT1 



Wakes of ram pressure stripped disc galaxies 7 




Figure 10. Velocity components in the wake, for medium ram- 
pressure, Mach number 2.53, i = 30°. Otherwise same as Fig. Isl 
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tively. For the projections along the 2-axis (right column) 
the different angles between the lower and upper line could 
be explained by the disc's rotation. In these panels, the disc 
is rotating counter-clockwise. Therefore, at the bottom side 
of the panel, the disc is rotating in wind direction, whereas 
on the upper side of the panel the disc is rotating antipar- 
allel to the wind direction. This leads to a slightly stronger 
stripping at the bottom side and a slightly stronger stretch- 
ing of the bottom side of the tail. For the projection along 
the a;-axis, the reason must be a different one because here 
the rotation does not induce a similar asymmetry. Here, the 
inclination makes the galaxy asymmetrical for the ICM flow. 
Moreover, the stripping proces s itself proceeds asymmetr i- 
cally for inclined galaxies (see iRoedieer fc Briiggeiil l2006ll . 
However, it is surprising that the asymmetry of the tail is 
similar for both inclinations and both projection directions. 

The fact that the widening of the tail proceeds in a very 
similar fashion independent of galactic inclination suggests 
that the dynamics of the widening are independent of the 
galaxy but intrinsic to the ICM flow. This also excludes the 
widening due to the rotation of the disc gas as presumed in 
the introduction as a major process. 

In the supersonic, medium ram-pressure case the tail 
widens, too, despite its irregular behaviour. The least widen- 
ing is observed in the supersonic high ram-pressure case. 
However, the velocity components perpendicular to the wind 
direction are high in the supersonic case, so that some widen- 
ing could be observed if the tail is traced to even larger 
distances behind the galaxy. 



Figure 11. Velocity components in the wake, for strong ram- 
pressure, Mach number 0.8, i = 30°. Otherwise same as Fig. Isl 



4 DISCUSSION 



Figs. |3| and |1] envelop the tail. They were fitted by eye. In 
each figure, they are identical for all panels, i.e. they are 
independent of time and projection direction. 

The difference in the angles of the envelopes between 
the two inclinations is very small. However, for the more 
highly inclined galaxy, the enveloping lines are closer to- 
gether. For higher inclinations, the width of the tail at a 
given distance to the galaxy is smaller. In both inclinations, 
the upper envelope has an angle of 13° with respect to the 
y-axis. The angles between the lower envelope and the y-axis 
are 9° and 10.5° for the inclinations of 75° and 30°, respec- 




Figure 12. Velocity components in the wake, for strong ram- 
pressure, Mach number 2.53, i = 30°. Otherwise same as Fig. Isl 



4.1 Observable galactic tails 

4.1.1 General 

Given that the stripped gas remains in the HI temperature 
range, our simulations suggest that deep HI observations 
could observe long galactic tails. For galaxies that experi- 
ence a strong ram pressure and lose nearly their complete 
gas disc in a few 100 Myr, the tails could be detected with 
less sensitive observations. In such cases, the tail can even 
separate from the source galaxy (see bottom panels of Fig.|SJ 
and appear as a large isolated gas cloud. In addition to the 
spatial separation, the "tail" gas cloud will also be separated 
from its source galaxy in velocity space. This makes it very 
difficult to recognise the original connection between the gas 
cloud and the source galaxy. An HI survey of the Virgo clus- 
ter IpDavics ct al. 2004) has indeed found two candidates of 
HI clouds without optical counterparts. 

Our simulations predict tail widths in the range of 60 
to 100 kpc at non-projected distances of ~ 100 kpc behind 
the galaxy. If a galaxy moves subsonically, if it is exposed to 
a medium ram pressure and seen along a favourable line-of- 
sight, even regular structures similar to von-Karman oscilla- 
tions could be detected. However, some additional processes 
could modify the appearance of the tails: 

• For galaxies passing through cluster centres, the ICM 
wind will be variable. The decrease of ram pressure after the 
peak ram pressure and back-falling material could change, 
both, the tail structure and the structure of the velocity 
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component in wind direction. The velocity components per- 
pendicular to the wind direction should not be sensitive to 
this process. However, in observations the components into 
and perpendicular to wind direction may be difficult to dis- 
entangle. 

• We neglected cooling and thermal conduction. These 
processes are important for the question whether the 
stripped gas stays observable as HI. If a tail is observed 
during heating or cooling, the structure observed in HI 
may be completely decoupled from the overall structure 
of all stripped gas because some parts drop out of the 
relevant temperature range. Ho wever, the observation of 
lOosterloo fc van Gorkoiiil (120051) suggests that at least a 
part of the stripped gas can survive as HI for a few 100 Myr. 
More detailed studies are needed to investigate the role of 
thermal conduction. 

• We also neglected viscosity in our simulations. Due to 
numerical viscosity, our simulations match the case of mod- 
erate to high Reynolds numbers. Magnetic fields in the ICM 
may suppress viscosity significantly, which could lead to 
higher turbulent velocities and greater widening of the tail. 
On the other hand, a substantial amount of viscosity would 
smooth the turbulent structures that develop in our simula- 
tions, and galactic tails may be narrower than we predict. 

• Local ICM motions as found in many cluster simula- 
tions (e.g. iFaltenbacher et all 120051: iTakizawal l2005l) could 
deform and bend the galactic tails. 

4.1.2 The case of NGC 4388 

So far, the only example of a long tail of cool ram pressure 
stripped gas is the cas e of NGC 4388 in the Virgo Clus- 
ter (.Oostcrl oo fc van Corkom .2005') . This galaxy's tail has 
a length of ~ 120 kpc in the plane of the sky, the tail width 
fiares from about 15 kpc near the galaxy to about 40 kpc 
at a projected distance of 100 kpc behind the galaxy. The 
ridge of maximum fiux density suggests a slight S-shape of 
the tail. 

In order to disentangle the true properties of the wake, 
the velocity component in the plane of the sky would have 
to be known. The length of the tail indicates that this 
component is substantial. The radial velocity of NGC 4388 
with respect to the cluster mean is ~ 1400 kms"^ (see 
IVollmer fc Huchtmeierl [20031 and references therein). Even 
if this galaxy had no velocity component in the plane of the 
sky, its radial velocity suggests already a slight supersonic 
motion through the ICM, which would be enough to qualify 
this galaxy's velocity as typical for cluster galaxies. Conse- 
quently, a substantial velocity component in the plane of the 
sky as indicated by the long tail suggests that NGC 4388's 
true velocity with respect to the cluster mean is rather high. 
Thus, it is very likely that this galaxy is moving superson- 
ically through the ICM. Therefore, our simulations suggest 
that the implied S-shape in the tail is not due to regular 
von-Karman oscillations. The observed S-shape could also 
be coincidental and not part of a regular pattern. In our 
supersonic simulations situations like this occur, the super- 
sonic tails are far from being perfectly straight. NGC 4388's 
tail shape may also be caused by other processes discussed 
above. 

Even if the velocity component in the plane of the sky 
is comparable to the radial component (in the cluster rest 



frame), this implies that the true length of the tail is at least 
1.4 times as long as the projected one. If the velocity compo- 
nent in the plane of the sky is only ~ 500 km s ~ ^ , the true tail 
is about 3 times as long as the projected one. A reasonable 
assumption for the ICM density at NGC 4388's position is a 
few 10~^ cm~^. With an assumed velocity of ~ 1800 kms~^ 
it is likely that this galaxy experienced a ram pressure in 
the range we studied here. In our simulations, we can re- 
produce such very long tails. The details of the tail length 
may also depend on the exact stripping history, i.e. the time 
dependence of the ram pressure during the cluster passage. 
Thus, comparisons with simulations of realistic cluster pas- 
sages that also trace the evolution of the galactic tail may 
be used to infer the galaxy's proper velocity. 

The width of the tail does not depend on the projec- 
tion direction. Thus, the tail of NGC 4388 is narrower than 
the tails in our simulations. This difference may be caused 
by reasons discussed above. E.g. if the ICM is more viscous 
than in our simulations, the tail will be narrower. Alterna- 
tively, some parts of the tail may already have dropped out 
of the relevant temperature range in order to be observable 
in HI. Comparisons with simulations including these effects 
could give important insights on the degree of viscosity and 
thermal c onduction in the ICM. 

Oost erloo fc van GorkornI i2005ll have also measured 
the radial velocity in the tail of NGC 4388. Close to the 
galaxy, the radial velocity agrees with the galaxy's radial ve- 
locity. Then it decreases along the tail by about 550kms~^. 
This means that the stripped gas at the end of the tail has 
reached only about 40% of the galactic radial velocity with 
respect to the cluster mean. However, it takes some time un- 
til the stripped gas is decelerated sufficiently to fall behind 
the galaxy. In this point, our simulations agree with the ob- 
servations. Also in our supersonic simulations, the stripped 
gas does not reach the full ICM wind speed inside our sim- 
ulation box, and gas at true distances of 100 kpc behind the 
galaxy has reached only about 50% of vicm- Again, details of 
the velocity structure depend on the exact stripping history. 

In addition to t he velocity decrease, 
lOosterloo fc van Gorkoinl i2005l) have measured the velocity 
width in the tail. They derive values of 100 to 200kms~^. 
This is surprisingly low as one would expect the stripped 
gas to preserve the rotation it had inside the galaxy, i.e. the 
width should be at least twice the rotation velocity. Our 
simulations show at least this width, and the supersonic 
simulations even show a higher width. This is in a ccord 
with the simulations of IVollmer fc Huchtmeierl |2003) that 
produce velocity widths of ~ 4 00kms~^. The reason for 
this d iscrepancy might be that lOosterloo fc van Gorkoinl 
(|2QQ^ have measured the velocity only along the ridge of 
highest flux density. Thereby, they may have incidentally 
selected a special fraction of the tail, like gas that originated 
from a certain part of the galaxy only. 

4.1.3 Other examples 

To our knowledge, no other long galactic tails have been 
observed in HI. Some observations exist in other wave- 
lengths. E.g. for CGCG 97-073 and 97-079 in the cluster 
A13 67 Ha and ra dio continuum tails have been detected 
fsee lGavazzi et al]ri995 - 2001). Both galaxies show tails of 
~ 75 kpc length. The Hq tails are rather thin and straight 
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and show little substructure. The radio continuum tails show 
a some widening. It will be interesting to observe these 
galaxies and others in deep HI for further analysis. 



4.2 Metal distribution in the ICM due to ram 
pressure stripping 

Gas stripped from galaxies is an important source of metal s 
for the ICM JSchindler et alJl2005l : Ibomainko et alJl2005l) . 

Therefore, the evolution of galactic tails plays an important 
role in the distribution of metals in the ICM. In our sim- 
ulations, we find that the width of the tails behaves like 
w/d ^ 0.4 (see e.g. Fig.j^l, where w is the tail width and d 
the distance from the galaxy. However, our simulations trace 
only the first 260 kpc behind the galaxy and the width of the 
wake is not going to increase linearly with the distance from 
the galaxy. In a simple model, we can model the turbulent 
widening of the wake as a diffusion process with a diffusion 
constant given by D ~ ^al, where cr is a typical turbulence 
velocity and I a typical size of turbulence eddies. Our sim- 
ulations suggest velocity dispersions of a ~ 250 km s^^ for 
galactic velocities of ~ 800 km s"'^ to cr 500kms~^ for 
galactic velocities of ~ 2500 km s~^. The typical eddy size, 
I, in our simulations is ~ 50 kpc. Thus, the time to reach a 
certain width w is given by 



1 

11) 
600 Myr 



3 

I'll 



100 kpc y V250kms 



I 



50 kpc 



(2) 



During this time, the tail length increases by ~ vicmt. Con- 
sequently, the tail length d and its width relate as 
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(3) 



This estimate neglects that the turbulence decays as the 
gas streams away from the galaxy. The decay of turbulence 
makes, both, a and / smaller with increasing distance to 
the galaxy and therefore, the above estimate gives an upper 
limit on w. 

These estimates could be included in future simulations 
of metal enrichment of the ICM by ram pressure stripping. 
Some of these simulations inject metals provided certain 
analytical criteria for ram-pressure stripping are fulfilled. 
These simulations could be refined by imparting the stripped 
material with a velocity dispersion to produce more realistic 
tails. 



4.3 Heating of the ICM by turbulence in the 
wakes 

The ongoing search for a heat source in galaxy clusters that 
prevents catastrophic cooling in the cool cores makes it inter- 
esting to estimate the extent to which the motion of galaxies 
can heat the ICM. 

The kinetic energy of a galaxy can heat its ambient 



medium via three mechanisms: (i) the flow around the 
galaxy can lead to turbulence in the wake of the galaxy as 
seen in our simulatio ns, (ii) dynarn i cal frictioii accel erates 
the surrounding gas JOstrikejIlQgl iKim et aLlliool . and 
(iii) in case of supersonic motion, the ambient gas is heated 
by a bow shock. For the first two mechanisms the induced 
motions need to be converted to heat by some dissipative 
process. 

As the galaxy moves through the ICM, it produces a 
turbulent wake. The energy loss of the galaxy may be esti- 
mated by 



dE, 



■gal 



dt 



= — Pram^projl'ICM, 



(4) 



where A is the projected surface area of the galaxy. Using 
conservative values, i.e. Pram = 6.4x lO"^'^ ergcm"^, vjcm = 
1000 kms~^ and A = 47r(4kpc)^, we obtain a loss rate ~ 
10'*'^ ergs~^. This energy can be found as turbulent energy 
in the ICM. For about 100 galaxies within a cluster, this 
heating rate is comparable to the cooling rate by radiative 
losses. 

Also, this heating rate is consistent with what we find 
in our simulations. A rough inspection of the velocity dis- 
tribution, see Fig. |2| and 1121 shows that the ICM is stirred 
within a tube of diameter ~ 40 kpc with a velocity disper- 
sion of ~ 400 kms~^. The ICM has an average gas den- 
sity of 10~^* gcm~'^ and the galaxy moves with 1000 km s~^ 
through the ICM. This yields a rate of kinetic energy injec- 
tion of 10'*^ ergs~^, which is very similar to our estimate 
above. 

It is interesting to compare this result to the energy 
input by dynamical friction. However, it should be noted 
that our simulations do not include the effect of dynamical 
friction (DF) since self-gravity is not included. 



dE^, 



dt 



= ^^-'^ dt 



4-K lnAG>icMMg\i 



(5) 



where In A is the logarithm of the ratio of maximum and 
minimum impact parameters and Mgai is the mass of the 
galaxy. Using the same values as above, we obtain an energy 
loss of ~ 6 X lO*'^ ergs~^, which is somewhat above the 
values derived for ram pressure heating. However, it is not 
clear how efficiently the energy is converted into heat. The 
accelerated gas may lead to propagating waves in the ICM 
that still need to be dissipated. In contrast, in the wake 
behind the galaxy, the ICM develops small-scale turbulence, 
which will dissipate locally. 

Moreover, ram pressure heating scales with ^proj, 
whereas dynamical friction heating scales with M^^i. Hence, 
for the large number of small galaxies ram pressure heating 
may prevail. 

Our simulations indicate that a supersonic galaxy stirs 
the ICM much more than a subsonically moving galaxy. 
Hence, gas that cools significantly below the local galaxy 
velocity dispersion will be heated more efficiently. 

Compared to ICM temperatures, the ISM lost from 
galaxies by ram pressure stripping is cold lO'' K) and 
thus might be considered an energy sink for the ICM. How- 
ever, the energy that becomes available by decelerating this 
stripped gas can be used for its heating. Given that cluster 
galaxies move with typical velocities comparable to the ICM 
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sound speed, this energy approximately balances the energy 
needed to heat the stripped ISM to ICM temperatures. 



5 SUMMARY 

We ran high-resolution 3D hydrodynamical simulations of 
ram pressure stripping of a massive disc galaxy focussing 
on the evolution of the galaxy's wake. We combined the 
following parameters: a medium and a high ram pressure, 
one subsonic and one supersonic ICM velocity, inclinations 
of 30° and 75°. We traced the tail to a distance of 260 kpc 
behind the galaxy. 

In our simulations, the wakes reached typical widths of 
~ 60 to 100 kpc at distances of 100 kpc behind the galaxy. 
The widening proceeds systematically with opening angles 
< 23° in subsonic runs. Depending on the projection angle, 
subsonic wakes can show oscillations similar to von Karman 
vortex streets. The wakes in the supersonic cases also show 
widening, but in an irregular fashion. 

We compared our simulations with the observed tail of 
NGC 4388. We could explain some features, but our simu- 
lated tails are wider than this observed one. We discussed 
possible reasons for this difference. 

Finally, we provided suggestions how our simulation re- 
sults could be used in cluster simulations that aim at pre- 
dicting the evolution of the metal distribution in the ICM 
due to ram pressure stripping. 
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